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Abstract 

Supercritical accretion flows inevitably produce radiation-pressure driven out- 
flows, which will Compton up-scatter soft photons from the underlying accretion 
flow, thereby making hard emission. We perform two dimensional radiation hydro- 
dynamic simulations of supercritical accretion flows and outflows, incorporating such 
Compton scattering effects, and demonstrate that there appears a new hard spectral 
state at higher photon luminosities than that of the slim-disk state. In this state, 
as the photon luminosity increases, the photon index decreases and the fraction of 
the hard emission increases. The Compton y-parameter is of the order of unity (and 
thus the photon index will be ~ 2) when the apparent photon luminosity is ~30Le 
(with Le being the Eddington luminosity) for nearly face-on sources. This explains 
the observed spectral hardening of the ULX NGC1313 X-2 in its brightening phase 
and thus supports the model of supercritical accretion onto stellar mass black holes 
in this ULX. 

Key words: accretion, accretion disks — black hole physics — hydrodynamics 
— radiative transfer 
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1. Introduction 

It is long known that astrophysical black holes (especially, black hole binaries) exhibit 
several distinct spectral states. These are, in the order of increasing photon luminosities (see, 
e.g., Esin et al. 1997): a low/hard state with a hard, power-law-like emission component, a 
high/soft state with a soft, thermal emission component (e.g.,Tananbaum et al. 1972; Dolan 
et al. 1979), and a very high state showing both spectral components (e.g., Miyamoto et al. 
1991; Kubota and Done 2004). In addition, there seems to exist another soft spectral state 
at even higher photon luminosities, comparable to the Eddington luminosity (Lg), called the 
slim disk state (Abramowicz et al. 1988; Watarai et al. 2000, see Kato et al. 2008 for a review 
and references therein). Although recent observations suggest its existence (Vierdayanti et al. 
2006; Okajima et al. 2006), a robust observational proof of this state has not yet been obtained. 

In this respect it is interesting to note that Ultraluminous X-ray sources (ULXs), which 
are recently found successively in the off-center region of nearby external galaxies, also show 
both the soft thermal and the hard power-law spectral states (Makishima et al. 2000; Cropper 
et al. 2004). Some of them, such as IC342 X-l and X-2, even show spectral transitions between 
these two states (Kubota and Makishima 2001). Notably, the typical photon luminosities of 
ULXs range between L p h ~ io 39-41 erg s -1 , which exceeds the Eddington luminosity for neutron 
stars and stellar-mass black holes. There are two possible models considered to account for such 
large photon luminosities: subcritical accretion (i.e., accretion below the Eddington accretion 
rate) onto an intermediate-mass black hole (IMBH, Miller et al. 2003; Miller et al. 2004; Cropper 
et al. 2004; Berghea et al. 2008) and supercritical accretion onto a stellar-mass black hole (King 
et al. 2001; Ebisawa et al. 2003; Okajima et al. 2006; Tsunoda et al. 2006; Vierdayanti et al. 
2006; Poutanen et al. 2007; Vierdayanti et al. 2008). Since the black hole masses of ULXs 
are poorly known, we cannot discriminate between these two models at present. The study of 
spectral states and spectral transitions, in analogy to black hole binaries, are one possibility to 
resolve this issue, but the situation turns out to be very complex. If the two states of ULXs 
correspond to the low/hard state and the high/soft state, the photon luminosities of ULXs 
should be sub-Eddington and thus the IMBH hypothesis will be favored. If, on the other hand, 
the soft state of ULXs corresponds to the slim-disk state, these ULXs should be supercritical 
sources, which means that the black holes cannot be IMBHs. The proper identification of the 
spectral states is thus essential for understanding the nature of accretion flows and black holes 
in ULXs. 

Interesting trends have been reported recently. Some ULXs (such as HoIX X-l) show 
spectral softening when the photon luminosity increases (Kajava & Poutanen. 2008). It is 
thus natural to conclude that the brighter, soft state corresponds to the slim-disk state. Then, 
its photon luminosity should be close to the Eddington luminosity, thereby supporting the 
supercritical hypothesis for ULXs. To our surprise, however, there exists another class of ULXs 
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(such as NGC1313 X-2) which exhibit the completely opposite behavior; that is, they show 
spectral hardening as the photon luminosity increase (Roberts et al. 2006; Mizuno et al. 2007). 
It is then reasonable to conclude that the observed soft-to-hard transitions correspond to the 
transition from the high/soft state to the very high state in the black hole binaries. In this 
case, the photon luminosities should be sub-Eddington and, hence, the IMBH model will be 
favored. However, we can propose another interpretation. If a hard spectral state exists in 
the supercritical regime and if these sources are in such a very bright phase, the supercritical 
model will survive. From a theoretical point of view, this seems to be a natural consequence: 
supercritical accretion flows inevitably produce radiation-pressure driven outflows and such 
outflows will Compton up-scatter soft photons, thereby making a hard emission component. 
The higher the photon luminosity is, the harder emission we expect. In this paper we will 
demonstrate that this is indeed feasible based on new two-dimensional radiation-hydrodynamic 
(RHD) simulations, which incorporate the Compton scattering effects. 

Two-dimensional RHD simulations of supercritical accretion flows around stellar-mass 
black holes were pioneered by Eggum et al. (1988) and followed by several authors (Okuda 
2002; Okuda et al. 2005; Ohsuga 2007a). Supercritical accretion in a quasi-steady state was 
first calculated by Ohsuga et al. (2005), who found that the total photon luminosities can 
indeed exceed three times the Eddington luminosity, whereas the apparent photon luminosities 
become more than ten times larger than the Eddington luminosity for a face-on observer (see 
also Heinzeller et al. 2006). They have also shown that hot outflows with their gas temperature 
higher than 10 9 K appear above the disk. This suggests that the inverse Compton scattering of 
soft disk emission by such high temperature plasmas should be important, though it was not 
yet studied quantitatively. This motivated us to perform extended simulations that incorporate 
Compton scattering effects. We find that the higher the mass accretion rate is, more high- 
temperature gas exists as outflow and, therefore, the larger becomes the Compton ^-parameter 
(see equation (6) in section 3 for the definition of the Compton y-parameter) . We thus expect 
that SED becomes harder when the photon luminosity is higher. This new state will explain 
the SED variation of some ULXs. 

The plan of this paper is as follows: In §2 we present basic equations and numerical 
methods. In §3 we present the structure of supercritical accretion flows with Comptonizing 
outflows and compute the Compton y-parameter. Finally, §4 is devoted to discussion of our 
findings. 

2. Methods 

Following Ohsuga et al. (2005), but also considering the effects of Compton scattering in 
the energy exchange between photons and electrons, we solve the RHD equations in spherical 
coordinates (r, 9,<f>). The radiative transfer equation is solved using the flux-limited diffusion 
approximation (Levermore & Pomraning 1981; Turner & Stone 2001). The general relativistic 
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effects are incorporated by a pseudo- Newtonian potential, ^ = —GMj (r — r s ) (Paczyrisky & 
Wiita 1980), where r s (= 2GM/c 2 ) is the Schwarzschild radius, G is the gravitational constant, 
M is the mass of the black hole (M = 10M Q is employed), and c is the speed of light. We 
assume that the flow is non-self gravitating, axisymmetric with respect to the rotation axis 
(i.e., d/d(j) = 0), and symmetric relative to the equatorial plane (where 9 = tt/2). We also adopt 
the a viscosity prescription (Shakura & Sunyaev 1973) and set a = 0.1. 

The basic equations are the same as those in Ohsuga et al. (2005) except that we add 
Compton heating/ cooling terms in the energy equations of gas and radiation; 

— + V-(eu) = -pV-v - AixkB + ckE + $ vis - r Co mp , (1) 

OE 

+ V-(E v) = -V-F - Vu:P + AtikB - ckE + r Comp . (2) 

Here v = (v r , v$, v$) is the velocity, p is the gas pressure, e is the internal energy density of 
the gas, B is the blackbody intensity, E is the radiation energy density, where the suffix 
represents quantities in the comoving frame. F = (F£, Ffj) is the radiative flux, P is the 
radiation pressure tensor, r Comp is the energy transport rate from the gas to the radiation field 
via the Comptonization, k is the absorption opacity, and $ v j S is the viscous dissipation function. 
We write r Comp as 

r Comp = 4w MT ~- T " d) ( E , (3) 

m e c z \ m pJ 

where T gas is the gas temperature, T ra d [= {Eq/o) 1 ^} is the radiation temperature, a is the 
radiation constant, p is the mass density of the gas, <tt is the cross section for electron scattering, 
m e is the electron mass, and m p is the proton mass. We assume a one-temperature plasma 
in which the electron temperature equals the ion temperature. This equation is obtained 
by integrating the Kompaneets equation over the frequency (see e.g., Padmanabhan 2000 for 
detail). 

Our methods are the same as Ohsuga (2006) except that we solve the energy equations 
including the Compton heating/cooling terms. We employ the operator-splitting method in 
the energy equations and additionally solve the following equations; 

Oc 

7^ Pcomp i (4) 

dE °-r (K\ 

-gj- — J- Comp , \0) 

using Newton-Raphson iteration, with bisection when Newton-Raphson method fails. 

We start the calculations with a hot, rarefied, optically thin atmosphere (i.e., without 
an initial optically thick disk). We solve the RHD equations numerically by an explicit-implicit 
finite-difference scheme on an Eulerian grid. The computational domain size is 3r s < r < 500r s 
and < 9 < ir/2. We imposed an absorbing boundary condition at r = r in = 3r s by attaching 
a damping layer, in which the physical quantities q gradually approach the inital values go as 
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<?new(i — n) = q(i) — (q(i) — qo)f(n), where r(i) = r in . The function /(n) is a smooth function 
which monotonically increases from to 1 as n increases. We continuously add mass through 
the outer boundary (r = r out = 500r s ) near the equatorial plane (0.457T < 9 < 0.5n) at a 
constant rate mi npu t. Here, rhi nput is the mass input rate normalized by the critical accretion 
rate, M ait =L^j c 2 . The injected matter is assumed to have a specific angular momentum 
corresponding to the Keplerian angular momentum at r = 100r s . On the other hand, we allow 
matter to escape freely but not to enter the computational domain at r = r out and < 9 < 0.457T. 
The number of grid points is (N r , Ng) = (96, 96). The grid points in the radial direction are 
distributed such that Aln r = constant, while the grids in 9 direction are distributed in such a 
way that Acos^ = 1/Ng. 

3. Numerical Results 

3.1. Quasi-Steady Structure 

We first overview the accretion flow and outflow properties. The overall evolution can 
be classified into two phases: a transient accumulation phase and a subsequent quasi-steady 
phase. In the former one, matter injected through the outer disk boundary creates continuous 
gas inflow and accumulates around r = 100r s , where the angular momentum of the injected gas 
equals that of the Keplerian rotation. Viscous processes allow the angular momentum of the 
gas to be transported outward, which drives the gas inflow. Eventually, the gas falls onto the 
black hole in a quasi-steady fashion. In this subsection, we fix the mass input rate at the outer 
boundary to a value rh input = 10 3 . 

Figure 1 shows the color contours of the gas temperature (top panels) and of the mass 
density overlaid with the fluid velocity vectors (bottom panels) in the quasi-steady state for 
models with and without Comptonization in the left and right panels, respectively. We easily 
notice that the gas temperature of the outflow, which was originally 10 9 -10 n K (Fig. lb), is 
now reduced considerably to ~ 10 75 -10 8 K, when we take into account the Comptonization 
effects (Fig. la). The gas temperature in the disk region, on the other hand, does not change 
appreciably. This is because the radiation temperature and the gas temperature are nearly 
equal due to frequent absorption and emission of photons in the disk. The ratio of the Compton 
cooling timescale (tcomp~e/|rc m P |) to the cooling timescale by the free-free and bound- free 
emission (t ffjbf ~e/|47rK£> — ckE \) is £com P Aff,bf~10 _5 in the outflow region and ~1 in the disk 
region. We also note that in the outflow region the ratio of the Compton cooling timescale to 
the escaping timescale of the outflow (t csc ~r out /0.1 c) is ^Comp/^csc^l0 ^ although that of the 
cooling timescale by free- free and bound-free emission to the escaping timescale is ts t bt/t esc ~10 1 . 
Therefore the outflowing gas is efficiently cooled by the Comptonization before the gas escape 
from computational domain (by contraries, the cooling by free-free and bound-free emission is 
not effective in the outflow region). 
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The outflowing gas is less dense in models including the Comptonization. A reduction 
in the outflow rate leads to an increase of the mass accretion rate onto the black hole and leads 
to a geometrically thicker accretion flow. 

We find that the accretion rate onto the black hole and the total viscous heating rate in 
the model with Compton cooling are larger than those in the model without Compton cooling 
by about 70% and 20%, respectively. The disk scale height is decreased by the Compton 
cooling in the outflow region. Therefore, the gas above the disk accumulates towards the disk 
plane (The radiative force in the disk region is reduced due to photon diffusion processes). 
Compared with the case without Comptonization, the amount of the radiation-pressure driven 
outflow is smaller, the mass accretion rate is larger, and the viscous heating rate is also larger 
in the case with Comptonization. The photon luminosity remains about the same, however. 
The energy of the photons which are swallowed by the black hole per unit time in the model 
with Comptonization is three times as large as that without Comptonization. The increase in 
the mass accretion rate makes photon-trapping effects more significant (Ohsuga et al. 2002). 
Most of the photons which are produced additionally by the inclusion of the Compton effects 
are generated around the innermost region of the accretion disk (i.e., in the photon-trapping 
region). Thus, they are swallowed by the black hole without escaping to the outflow region 
(Ohsuga and Mineshige 2007). This is the reason why the observed photon luminosity does not 
change appreciably by the inclusion of the Compton effects. 

3.2. Hardness of the Photon Spectrum 

The Compton cooling of the outflow region affects the radiation spectra. In this subsec- 
tion, we use time-averaged simulation results to discuss the spectral properties in terms of the 
Compton y-parameter, which indicates how much seed soft photons are up-scattered by the 
hot electrons in the outflow. The y-parameter is given by 

y = 2 max(r es , r es ) , (6) 

m e cr 

where r es is the Thomson optical depth of the outflow. Note that we assume a one-temperature 
plasma and, hence, electron and proton temperatures are the same in this paper. 

We calculate the Compton y-parameter in the following way. First, we specify a ray, a 
straight line that connects a point at the outer boundary (r ont ,8) and the origin (r = 0) with an 
inclination angle of 9 (measured with respect to the z-axis). Second, we define a photosphere 
for that ray by the position where t c r{6) = \Jrj9)(jj9) + r cs (9)) = 1. Here the optical depth 
for absorption r a (#) and for electron scattering T es (8) are integrated from the outer boundary 
of the computational domain. This segment between the photosphere and the outer boundary 
is used to compute the ^/-parameter. If the medium along the ray is optically thin, we integrate 
the opacity inwards to the inner boundary. Third, we obtain the opacity-weighted mean plasma 
temperature along the ray outside the photosphere by 
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^ = ET gas (n,9)Ar es (n,9) 

where Ar es (ri,9) is the electron scattering opacity in the i-th grid point along the ray, and 
r cs = 12A-T cs (r[,8). By separating the summation along the ray to the high temperature region 
up-scattering photons (T e > 10 7 K) and the low temperature region down-scattering photons 
(T c < 10 7 K), we confirm that the contribution to the y-parameter from the low temperature 
region is less than 1. Thus we include only the grid points where T gas (ri,0) > 10 7 K to compute 
T e , r es , and y in equation 6. Otherwise, y will be overestimated. Figure 2 shows the y-parameters 
in the outflow region. The solid curves show the results for models with Comptonization for 
various mass input rates, whereas the dashed curve shows the ^-parameter for a model without 
Comptonization. When we include Comptonization effects in the RHD simulation, the y- 
parameter decreases from ^10 3 to 10~ 2 -10° because the gas temperature and the gas density in 
the outflow region decrease (see Fig. 1). Importantly, we confirm that the Compton ^-parameter 
in models with Comptonization is consistent with the observations of ULXs, while that in 
models without Comptonization is too large to explain them. We also find that the ^/-parameter 
increases with an increase of rhmput (and, hence, an increase of the photon luminosity). This is 
because the larger m inpu t is, the larger becomes the electron number density, and, hence, the 
larger becomes the Thomson optical depth. The SED becomes harder as the photon luminosity 
increases. 

Table 1 summarizes the photon index T, the ^-parameter, the Thomson optical depth 
r cs , the averaged gas temperature T e , and the normalized isotropic (apparent) photon luminos- 
ity L p h(iso)/-^E, for a single ray with an angle of #/(7r/2)~0.05 (i.e., nearly face-on view). Here 
we calculate the isotropic photon luminosity by £ ph ( iso )=47T (r out ) 2 FQ (rout,^). The photon index 



T (=— l/2 + y9/4 + 4/?/) resulting from the unsaturated Comptonization (Rybicki & Lightman 
1979) is computed using the y-parameter at 6'/(7r/2)~0.05. In Table 1, we also present the 
normalized kinetic luminosity L^/L^, the normalized photon luminosity L^/Le, the mass ac- 
cretion rate onto the black hole m(=M/M cri t), the mass outflow rate m ou tflow(=^outflow/-^crit), 
which are integrated over the all angle. The kinetic luminosity is defined as the mechanical 
energy of outflow which exceed the escape velocity v csc at the outer boundary (r = r out ): 

L kin =(r out ) 2 J ^p(r out ,9)v r (r out ,9) 3 dtt. (8) 
The mass accretion rate is defined at the inner boundary as 

M=(r in ) 2 | p(r in ,#)max[-t; r (r in ,fl),0]da (9) 
The mass outflow rate is evaluated at the outer boundary: 

M outaow =(r out ) 2 J p(r out , 6) max[t; r (r out , #), 0]dlX (10) 

Figure 3 shows the dependence of the space-averaged gas temperature of the outflow, the 
Thomson optical depth, the ^-parameter, and the photon index for an observer located in the 
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direction of #/(7r/2)~0.05 on the isotropic photon luminosity. We find a negative correlation 
between T and L p h(iso)- This is because max(r es , r e 2 s ) increases with an increase of fa. 

4. Discussion 

As mentioned in § 1, we expect a hard spectral state exists in the supercritical regimes 
for the supercritical accretion model for ULXs (i.e., central black holes are stellar-mass black 
holes) to be viable. In this paper, we demonstrated that such a hard state can really appear 
in the case of the supercritical accretion onto black holes as a result of the combination of an 
optically thick, slim-disk type inflow and a Comptonizing outflow. In this new supercritical 
state, higher photon luminosities correspond to harder SEDs, in analogy to black hole binaries 
in the very high state, but at even higher photon luminosities. 

Figure 4 illustrates schematically our proposed scenario for the variations in the disk 
geometry and the spectral properties according to the changes in the mass input rate. From 
top to bottom: (a) Comptonizing outflow state, (b) slim disk state, (c) very high state, and (d) 
high/soft state. In the slim disk state, in which the photon luminosity barely exceeds Le, the 
amount of gas outflow is negligible. Hence, a pure thermal emission component is detectable. 
However, as the mass input rate increases, the amount of radiation-pressure driven, mildly 
hot (~ 10 7 ' 5_8 K) outflow increases, thereby more soft photons from the accretion flow being 
Compton up-scattered. The higher the accretion rate is, the harder becomes the SED. This 
transition can account for the spectral transition reported for some ULXs (e.g., NGC1313 X-2), 
which show a positive correlation between the SED hardness and the photon luminosity. 

What will be an observational test of our proposed spectral state? Unfortunately, line 
absorption by the outflow material is not expected, since the temperature of the Comptonizing 
outflow is so high that all the metals in the Comptonizing outflow are completely ionized. 
Instead, we might search for sources which show both slim-disk features and a Comptonized 
blackbody component. The slim-disk features can be evaluated through the spectral fitting 
by using the extended disk-blackbody (DBB) model (Mineshige et al. 1994), in which the 
temperature profile is assumed to be proportional to r~ p , with r being the distance to the 
central black hole and p being a fitting parameter. The p-value depends on the disk model: 
p=0.75 in the standard-type disk and p=0.5 in the slim disk (see Chap. 10 of Kato et al. 2008). 
Using this technique, Vierdayanti et al. (2006) found that at least some ULXs are powered by 
supercritical accretion onto stellar-mass black holes. If there are ULXs which show both a soft 
component with p = 0.5 and a Comptonized hard component, they will provide strong support 
for our model. 

Next, let us discuss the effect of the bulk (motion) Comptonization. The bulk 
Comptonization becomes significant when the speed of the bulk motion becomes rela- 
tivistic. The bulk Comptonization dominates the thermal Comptonization when £ = 
v/tUkBT^/rric) 1 / 2 > 1 (Blandford and Payne 1981). In our simulation results, the outflow 
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velocity is ~ 0.5 c, while the gas temperature is ~ 10 7 - 5 " 8 K, giving rise to £ ~ 1 in the region 
near the polar axis. This means that the bulk motion Comptonization is not negligible, com- 
pared with the thermal Comptonization. We calculate the ^-parameter of bulk Comptonization 
by Vbuik = {m e v 2 /3m e c 2 )max(r es , r e 2 s ), finding that bulk Comptonization becomes comparable to 
thermal Comptonization in the region around the polar axis. Hence, we expect substantial bulk 
Comptonization effects in the spectra of nearly face-on sources. Caution should be taken here, 
however. Our simulations are likely to overestimate the bulk velocity of the outflow, since the 
effect of the radiation drag, which decelerates the hydrodynamic motion, is not considered in 
our code. On the other hand, we evaluate the radiation temperature as T ra d = (Eo/a) 1 ^ 4 , which 
results in underestimation of the radiation temperature, when the radiation spectra become 
harder than a Planck distribution. Therefore, we may overestimate the cooling rate of the gas 
via Comptonization, and, hence, underestimate the gas temperature. To evaluate the effects 
of bulk Comptonization more accurately, full non-gray RHD simulations are required. Future 
advances in computer technology will enable the study of full non-gray RHD simulations of 
supercritical accretion flows. We should also note that exponential 'photon breeding' mecha- 
nism, which converts a kinetic energy of a relativistic jet into a radiation energy via the bulk 
Comptonization by electron-positron pair plasmas in a jet (this is a positive feedback mecha- 
nism because pair plasma is created by high energy photons), is not efficient. This is because 
the outflow speed obtained in this study is about 0.5c although exponential photon breeding 
requires highly relativistic motion whose Lorentz factor T exceeds 3-4 (Stern & Poutanen 2006). 

We showed that in supercritical accretion flows , the kinetic luminosity of the outflow Lki n 
is 0.075, 0.30, 0.55, 1.1 L E for m input = 3xl0 2 , 10 3 , 3xl0 3 , 10 4 respectively (see Table 1) and 
that the momentum flux of the outflow is ~ 2L^\ n /v ~ 0.2 — 2 L^/v ~ Le/c, where v (~0.5c) is 
the outflow velocity. The mass, momentum, and energy outflow from the accreting black holes 
can affect the surrounding interstellar gas. This feedback has been discussed in supermassive 
black holes (AGN feedback) and should also be relevant in Galactic microquasars. By assuming 
that the momentum flux of the outflow is ~ L E /c as suggested by King and Pounds (2003), 
King (2003, 2005) showed that outflows from a supercritically accreting black hole at the center 
of the galaxy create a bubble, inside which stars are formed. This mechanism can explain the 
relation between the black hole mass and the velocity dispersion of the bulge [M — a relation). 
The results of our simulations for mi nput = 10 3 and 3xl0 3 are consistent with the assumptions 
of King (2003, 2005), if the black hole at the galactic center is accreting the gas supercritically. 
We note that compared to Ohsuga (2007b), which performed two-dimensional RHD simulations 
without Comptonization, higher mass input rate is required in our simulations to launch the 
outflow with the momentum flux ~ L-^/c. This is because the mass outflow rate decreases when 
Comptonization is included, as we mentioned in §3.1. 

In this paper we simply assumed the ^-prescription of the viscosity. However, it is now 
widely accepted that the MHD turbulence driven by the magneto-rotational instability (MRI) 
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provides a promising mechanism for the angular momentum transport in high-temperature 
flows (Balbus and Hawley 1991; Hawley et al. 1995; Balbus and Hawley 1998). In addition, the 
Parker instability, driven by the magnetic buoyancy (Parker 1966), creates magnetic loops in 
the disk corona (Machida et al. 2000). In order to investigate supercritical accretion phenomena 
without introducing the a- viscosity, we need to perform global radiation magnetohydrodynamic 
(RMHD) simulations (such simulations are pioneered by Ohsuga et al. (2009), but Compton 
effects are not included). To treat the relativistic effects of RHD such as radiation drag, it is 
also necessary to extend our numerical model of the radiative transfer to be valid up to v jc ~ 1. 
This remains as future work. 

We would like to thank T. Hanawa, S. Hirose, T. N. Kato, H. Oda, T. Sano, H. Takahashi, 
and S. Takeuchi for helpful conversations. Numerical simulations were carried out on Cray XT4 
at Center for Computational Astrophysics of National Astronomical Observatory of Japan and 
on Express5800/120Rg-l at Cybermedia Center Osaka University. This work is supported 
in part by Ministry of Education, Culture, Sports, Science, and Technology (MEXT) Young 
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with Comptonization 



without Comptonization 




R/r a Rfr B 

Fig. 1. Distribution of temperature at t = 40 s (top) and mass density averaged dur- 
ing t = 30-50 s (bottom). (Left) The model including Comptonization. (Right) 
The model without Comptonization. Arrows represent velocities of fluid motions. 
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Fig. 2. Compton y-paramctcr in the outflow region with respect to the polar angle 6/(tv/2) for 
"linput = 3xl0 2 (red), 10 3 (green), 3xl0 3 (blue), and 10 4 (magenta). Solid curves represent re- 
sults for models with Comptonization , the dashed curve shows the results for the model without 
Comptonization. The y-parameters become zero when 6/(ir/2) > 0.12 and > 0.22 for mi npu t = 3xl0 3 
and 10 4 , respectively. This is because the gas temperature of the outflow is below 10 7 K. 



Table 1. Mass accretion rate, mass outflow rate, luminosities and properties of Comptonizing outflows.* 



^input 


771 


''"outflow 


Lkin/ Le 




-^ph(iso)/ Le 


max(r os ,r 2 s ) 


T e 


y 


r 


3xl0 2 


1.3xl0 2 


3.7X10 1 


0.075 


2.4 


3.6 


1.5xl0~ 2 


3.7xl0 7 


3.8xl0" 4 


100 


10 3 


2.4xl0 2 


7.5xl0 2 


0.30 


3.2 


11 


3.9X10- 1 


4.2xl0 7 


l.lxlO" 2 


19 


3xl0 3 


4.9xl0 2 


2.5xl0 3 


0.55 


4.5 


21 


4.4x10° 


8.3xl0 7 


2.5xl0 _1 


3.8 


10 4 


1.3xl0 3 


8.6xl0 3 


1.1 


8.8 


39 


l.lxlO 2 


6.0xl0 7 


4.4x10° 


1.3 



* Minput is mass input rate, m is mass accretion rate, m ou tflow is mass outflow rate, Lkin is kinetic luminosity, 
L p h is photon luminosity, L p h(iso) is isotropic photon luminosity, Le is Eddington luminosity, r cs is Thomson 
optical depth, T c is electron temperature, y is Compton y-parameter, and T is photon index, respectively. We 
note that m, m outn ow7 Lkin, and L pn are integrated over the all angle. By contrast, L p h(iso) is evaluated by 
iph(iso) = 47r(r out ) 2 -FQ (r out ,0), where #/(7r/2)~0.05 The quantities max(r es ,r 2 s ), T e , y and F are also evaluated 
at 0/(tt/2)~O.O5. 



13 



I - 



a" 

1 




15 2H 25 

■^ixi/^E 

Fig. 3. The gas temperature of outflow averaged along the ray outside the photosphere T (top panel), 
y-parameter, Thomson optical depth max(T cs ,T c 2 s ) (middle panel), and the photon index T (bottom panel) 
as functions of photon luminosity Lp^/L^ at #/(7r/2)~0.05. (Cross: mi nput — 3xl0 2 , Circle: TOi npu t = 10 3 , 
Square: rhmput = 3xl0 3 , Triangle: 77ii, lpu t = 10 4 ). The dotted line in the middle panel represents y = 1. 
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Fig. 4. Schematic pictures of the states of accretion disks. From top to bottom: (a) 
Comptonizing outflow state, which consists of supercritical accretion flows (i.e., slim disks) 
and Comptonizing outflows, (b) slim disk state, (c) very high state, and (d) high/soft 
state, The SED of the new state is harder than that of the slim disk state, because the 
new state includes a hot outflow which up-scatters seed photons from the underlying disk. 
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